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ABSTRACT 

The intervening z — 0.7443 Mgll absorption system in the spectrum of MC 1331 + 170 shows an unusual series 
of line pairs, each with Av ~ 30 km s _1 . These lines could be explained as the shells of expanding superbubbles 
residing in the outer regions of an edge-on spiral galaxy visible in the optical image of the MC 1331 + 170 
field. The color and brightness of this galaxy make it the most likely candidate z = 0.7443 absorber, though two 
other galaxies in the quasar field could also be contributing to the Mgll absorption profile. Kinematic models of 
absorption from compact groups and galaxy pairs produce profiles largely inconsistent with the observed Mg II 
spectrum. Superbubbles would naturally generate more regular structures such as those observed. Photoionization 
models of the superbubble shell are consistent with the observed profile for many realistic physical conditions. 
In a pure superbubble model, the large velocity spread of the Mg II absorption system is inconsistent with the 
expected spread of a quiescent, rotating disk. This requires unusual kinematics within the host galaxy, perhaps 
due to a recent interaction. 

Subject headings: quasars: absorption lines — interstellar medium — superbubbles 



1. introduction 

Strong Mgll absorbers show complex kinematic structure 
consistent wit h a mixture of radial infa ll from the halo and a 
rotat ing disk (tLanzetta & Bowen 19921; Charlton & Churchill 
1998; phurchill & Vogt 2001 ). A large variety of phenomena 
can contribute to the profiles, such as Hi regions, Hll regions, 
high-velocity halo clouds, and superbubbles. Strong Mgll ab- 
sorb ers almost always have Lb > 0-08 L%, galaxies wit hin 35 
kpc dSteidel, Dickinson, & Persson 1994tpteidel 1995b - One 
particular absorption system, at z — 0.7443, in the spectrum 
of MC 1331 + 170 shows multiple pairs of lines, each with 
Av ~ 30 km s _1 . There are four strong optical sources within 
5" of the quasar, one of which is an edge-on spiral. These are 
the only detected optical sources which could have impact pa- 
rameters less than 35 kpc. 

OB associations are capable of producing large shells of 
swept up ISM using the energy input of stellar winds and 
multiple supernovae. Such structures, k nown as superbu bbles, 
have been obser ved in nearby galaxies (|Kamphuis 1993t Sun 



geon et al. 1999) as well as in our own galaxy ( [Heiles 1979 



McClure-Griffiths et al. 2000[ ). Mg II and Fe II absorption would 
arise within the cold superbubble shell, appearing as pairs of 
lines from the oppositely expanding sides of the shell intersect- 
ing the line of sight. 

In our galaxy, superbubbles are detected in Hi as filaments, 
with inferred radii ranging from 100 pc to 2 kpc (e.g. H eiles 
1979). Extragalactic superbubbles are manifest as holes in the 



Hi map of face-on galaxies or as Hi emission-line filaments 
in edge-on galaxies. Single structures of this type are observed 
with sizes similar to Galactic bubbles and masses between 10 7 
and 10 s Mq. Smaller shells (even single-star stellar wind bub- 
bles) are also thought to exist, but cannot be detected in ex- 
ternal galaxies. Superbubbles larger than 1 kpc in radius are 
thought to be formed by either a merging of separate shells or 
by a proce ss of self-propagating star formation (e.g. B omans 
et al. 1996). The number of bubbles in a given galaxy varies 
dramatically from almost none to several hundred, perhaps de- 
pending on the galaxy's history of interactions. 

Stellar wind bubbles for a single early-type star were first 
modeled by Weaver et al. (1977). Mac Low & McCray (1988) 
extended this concept to clusters of O and B stars, proposing 
Galactic superbubbles as the explanation for the Hi filaments 
observed by Heiles (1979). According to these models, a super- 
bubble begins as an oversized stellar wind bubble in the ISM 
around an OB association. The swept up material from this 
cavity quickly collapses into a cold, dense shell with a density 
greater than four times that of the ambient ISM. The shell be- 
gins expanding with the onset of the first supernovae in the OB 
association, which contribute their energy and metals to the su- 
perbubble. During this, the adiabatic phase, the supernovae can 
be approximated as a continuous energy source and the bubble 
radius increases as i 3 / 5 . The interior is sparse, highly ionized, 
and isobaric during this phase. The continuing energy input 
maintains the expansion at a rate of t 1 / 2 , in excess of the rate 
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expe cted for a pure snowplow solution ( Mac Low & McCray 
1988). 

In this paper, we explore the hypothesis that the line splittings 
in the z = 0.7443 absorber arise from expanding superbubbles 
in the vici nity of an observed edge-on spiral galaxy (C hurchill 
et al. 1995). In § we present the HIRES/Keck absorption 
profiles of detected low ionization transitions and an archival 
//5T/WFPC2 image of the quasar field. In § ||J we constrain 
the redshifts and luminosities of the four galaxies in the quasar 
field using the galaxies' colors and magnitudes. We consider 



of D LAs (100 — 200 K) occur when they lie in spiral galaxies 
(e.g. Kanekar et al. 200C). Even if we are looking through the 
disk of a spiral galaxy (as we suspect), the 2 1-cm data do not re- 
strict the nature of the absorber. However, the Mg II absorption 
arising from known DLAs has the c haracteristic of being to tally 
saturated across the entire profile (Churchill et al. 2000b) and 
typically exhibits a total velocity spread of 50 — 100 km s _1 . 
These features are not present in this system, suggesting that it 
is not a DLA. 



groups give rise 
5L we provide a 



the possibility that galaxy pairs and compact 
to the observed Mg II absorption in § [5]. In § I 
description of the photoionization models used in testing the 
superbubble hypothesis. We discuss our results and summarize 
in§! 



2. 
2.1. 



DATA AND ANALYSIS 



HIRES/Keck Spectra 



16.7 quasar 
Keck Telescope's 



2.2. HST/PC2 Image 

Images of the MC 1331 + 170 field obtained through the 
F702W (R) and F814 W (I) filters were retrieved from the 
archive (Program 5351, |^e Brun et al. 1997 ). Pipeline-reduced 
images for each filter were averaged and cleaned using the 
IRAF task GCOMBINE. The task COSMICRAYS was used to 
remove hot pixels. To increase the signal-to-noise ratio for pre- 
sentation purposes, the R and I band images were combined. 
For our analysis, we use the photometry as presented in Sec- 
tion 2, Table 8 of Le Brun et al. 1997. Two realizations of this 
ffi>r/WFPC2 image, with the quasar centered in the PC, are 
presented in Figure |[ The left-hand panel has a larger dynamic 
range and shows three resolved objects within 5" of the quasar, 
labeled in order of angular distance from the quasar. An addi- 
tional source at ~ 0.8" is apparent upon point-spread function 
(PSF ) subtraction of the quasar (Figure 16, of Le Brun et al. 
1997). In the right-hand panel, G5 is emphasized. 

3. THE MC 1331 + 170 FIELD 

3.1. Procedure 

The photometry of Lebrun et al. (1997) gives a limiting mag- 
nitude of 26 in R, which implies that unresolved sources at 
z = 0.7443 brighter than — 0.04 L* B should be seen in the 
image. Because strong Mgll absorbers nearly always have 
Lb > 0.08 L* B cou nterparts within 35 kpc (S teidel, Dickin- 
son, & Persson 1994; [Steidel 1995 ), we expect to detect the 
galaxy giving rise to the z = 0.7443 absorption. 

There are four galaxies which could be within 35 kpc of the 
quas ar in the HST/WFPC 2 image (Figure^, G2-G5, numbered 
as in Le Brun et al. 1997). Inspection of the sources shows an 
edge-on spir al morphology for G 5 and an elliptical shape for 
G3 (also see Le Brun et al. 1997). G4 shows some structure 
(possibly irregular) and the location of G2 in front of the quasar 
makes morphology difficult to determine. 

Unfortunately, redshifts for G2-G5 have not been obtained. 
We must, therefore, determine how likely it is that each candi- 
date galaxy is associated with the absorber at each of the red- 
shifts. This will allow us to assess which galaxy/galaxies could 
give rise to the z = 0.7443 absorption. We do this by com- 
paring the apparent magnitudes and colors of each galaxy to 
redshifted model spectral energy distributions (SED) for dif- 
ferent morphological types. However, we appreciate that there 
may not be a one-to-one correspondence between absorbers 
and galaxies. Two or more of the observed galaxies could give 
rise to absorpt ion and lie at the same redshift ( Bowen, Blades, 
& Pettini 1995). Additionally, DLA absorption can arise from 
relatively low luminosity and/or low sur face brightness galaxies 
which are not seen in the i mage (e.g., |Rao & Turnshek 2000 ; 
Bouche et al. 200C; also see Bowen, Tripp, & Jenkins 2000). 

,5 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by AURA, Inc., under contract to the NSF. 

4 The Mgll at z = 1.3284 is blended with the Fell from the higher redshift systems (Steidel & Sargent 1992) and a precise equivalent width measurement is not 
possible, but we estimate VF r (2796) = 0.8 A. 



The spectrum of the z = 2.084, V 
MC 1331 + 170 was obtain ed with the W. M 
HIRES/Keck spectrograph ( |Vogtetal. 1994| ) on the night of 19 
March 1993. The slit width was 0.861" (yielding R = 45, 000). 
A 2 x 2 pixel binning (in both the cross-dispersion and disper- 
sion directions) was used during read out. The resulting resolu- 
tion is R = 34, 000 (v = 8.8 km s _1 per resolution element). 
The wavelength coverage is 4263 to 6733 A, with gaps above 
5 100 A, because the CCD did not capture the free spectral range 
of the echelle orders above this wavelength. The spectrum was 
reduced in the usual manner using IRAF 3 v2.10. All wave- 
lengths are vacuum and heliocentric corrected. 

There are four Mgll systems in the spectrum of MC 1331 + 
170 at z = 0.7443, z = 1.3280, z = 1.776, and z = 
1.786. Due to its unusual structure, our focus is on the z = 
0.7443 system. The detected transitions (Fell A2587, A2600, 
Mgll AA2796,2803, Mgl A2853) for the z = 0.7443 system 
are shown in Figure |l| The profile fits (solid curves) were ob- 
tained using Voigt profile decomposition. This yields the num- 
ber of components, their line-of-sight velocities, column den- 
sities, and Doppler b parameters. We used MINFIT , a x 2 min- 
imization code of our own design ( [Churchill 1997 ). The ticks 
above the normalized continua in Figure [j] give component ve- 
locity centroids. Table |l] lists the resulting parameters from the 
Voigt profile fitting. 

The equivalent width of the z = 0.7443 system is 
W r (2796) = 1.81 A. An HST/STIS spectrum that covers the 
ClV AA1548, 1550 doublet is available in the archive (Program 
7271). Unfortunately, the signal to noise is ~ 1 at the posi- 
tion of the Civ absorption, yielding no useful equivalent width 
limit. The other systems are also strong 4 . Damped Lya ab- 
sorption (DLA), N(Ul) > 2 x 10 20 c m' 2 is observed in the 
z = 1.776 system ( phaffee et al. 198 80 . 

Rao and Turnshek (2000) found that 50% of Mgll absorbers 
with W U2796) > 0.5 A an d VF r (2600) > 0.5 A are DLAs 
(also see Boisse et al. 1998[ ). This criterion is satisfied by the 
z = 0.7443 system, but the Lya of this system is below the Ly- 
man limit break of the z = 1 . 776 DLA and cannot be measured. 
Neut ral hydroge n 2 1-cm measurements of the z = 0.7443 sys- 
tem ( |Lane 2OOC0 give iV(Hl) < 1.5xl0 18 (T s ) cm~ 2 , where 
(T s ) is the spin temperature. The minimum spin temperatures 
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In order to obtain the Lb values expected for each galaxy 
at each redshift, we first obtained B and R magnitudes from 
model SEDs for a set of standard morphologies (B ruzual & 
Chariot 1993; |Kinney et al. 1996[ ). Models for E, SO, Sab, 
Sbc, and Scd galaxies (/x = 0.95,0.30,0.20,0.10,0.01, re- 
spectively 5 ) were modelled with a 16 Gyr stellar population 
with exponentially decreasing star formation. A redder ellip- 
tical galaxy (E2) was modeled using a composite spectrum of 
observed ellipticals from Kinney et al. 1996. For an irregu- 
l ar galaxy (Im), the obse rved SED of NGC 4449 was adopted 
(Bruzual & Chariot 1993). Two starburst models were also con- 



sidered, both with a constant star formation rate, assuming a 
Salpeter IMF. The ages of the starburst galaxies (Starburst 1 and 
Starburst 2, Bruzual & Charlaut 1993) are 10 6 and 10 7 years. 

The K-corrections were performed by redshifting the model 
SEDs to the absorber redshifts and convolving them with the 
F702W and F814W filter response curves. A zero-point nor- 
malization was performed using the measured R band appar- 
ent magnit udes listed in Table H (also see Table 8 in L e Brun 
et al. 1997). Finally, the Mb values were calculated using 
Ho = 50 km s _1 Mpc -1 and qo = 0.05. Figure [}] shows 
R — I color vs. redshift for the standard morphological types. 
The dotted lines on the figure correspond to the redshifts of 
the Mgll absorptions systems in the MC 1331 + 170 spectrum 
and the dashed lines to the colors of the four galaxies, G2-G5. 
Intersections of the lines represent permitted locations on the 
diagram of the galaxies in the field and allow us to constrain 
their type at each possible redshift. 

3.2. Results 

Table || summarizes the permitted morphological types and 
their calculated Lb (in units of L* B , using M B = —20.9) at 
each absorber redshift (excepting z = 1.786 because of its 
proximity to z — 1.776). 

The errors on the color of G2 and magnitude might be quite 
large due to the PSF subtraction of the quasar. The color of 
G2 makes it unlikely to be at z — 0.7443, and it must be an 
elliptical if it is at z = 1.3284 or z ~ 1.8. If only its color 
is considered, it could be an E or SO at z ~ 1.8, but the de- 
rived Mb values give Lb ~ 15 L* B for an E and Lb ~ 5 L* B 
for an SO. Such luminous galaxies are exceedingly rare and we 
conclude that G2 is probably at z = 1.3284. 

The morphology and color of G3 seem to indicate that it is an 
early-type galaxy, but its redshift is difficult to constrain. If at 
z = 0.7443, G3 would have Lb ~ 0.1 L* B and could be an E, 
SO, or early-type spiral. Due to the low luminosity, it is uncer- 
tain whether G3 would be giving rise to strong Mg II absorption 
at this redshift. G3 could be an SO at either z = 1.3284 with 
L B ~ 1 L* B or z — 1.8 with L B ~ 4.3 L* B . If it is the latter, 
the galaxy would be unusually luminous. Therefore, it is most 
likely an E or SO at z = 0.7443 or z = 1.3284. 

G4 is very blue and would have to be actively star-forming 
to be at z = 0.7443. Just as with G3, however, it may not be 
luminous enough (Lb ~ 0.1 L* B ) at this redshift to be expected 
to give rise to strong Mg II absorption. The irregular structure 
observed in the image may indicate an Im galaxy. The color 
and luminosity (Lb ~ 0.1 L* B ) are also consistent with such 
a morphology at either z = 1.3284 or z ~ 1.8. The color- 
redshift curves converge toward R — I ~ at high redshift for 
the spiral types, allowing an Scd or Sbc identification for G4 at 
z-1.8. 

5 fi is the fraction of the mass in the stars after one Gyr (Bruzual 1983) 



The most striking galaxy in the MC 1331 + 170 field is G5, an 
edge-on spiral. This is our favored candidate for a z = 0.7443 
absorber, with a luminosity of Lb ~ 1 L* B . The low luminosi- 
ties of G3 and G4 make them marginal cases for producing such 
strong Mg II absorption. Additionally, the luminosity of G5 be- 
comes larger for the higher absorber redshifts (^2 — 3 L* B ); 
thus, it ~ 10 times less likely to be at the higher redshifts. The 
color of G5, however, does not constrain its redshift and allows 
it to be a star-forming galaxy at any of the higher redshifts. If 
this is so, G3 or G4 would be responsible for the z — 0.7443 
system, but their small luminosities at this redshift make this 
unlikely. 

In conclusion, G5 is the galaxy most likely to be producing 
z = 0.7443 absorption, but G3 and G4 could also be contribut- 
ing to the absorption at that redshift. If so, they would have 
small luminosities that make them border-line candidates for 
strong Mg II absorption. 

4. KINEMATIC MODEL SIMULATIONS OF GROUPS AND PAIRS 

Though G5 is the best candidate for the z — 0.7443 ab- 
sorber, G3 and G4 also have colors and luminosities consistent 
with this redshift. Furthermore, the blue colors of G4 and G5 
seem to indicate recent star formation, which could easily have 
been triggered by an interaction between two or more galaxies. 
As such, two compelling possibilities are that the line of sight 
passes through several members of a compact group of galax- 
ies or through a galaxy pair, giving rise to the observed Mg II 
absorption profile. 

Charlton and Churchill (1998) found that the kinematics of 
strong Mg II absorbers can be described by assuming that the 
absorbing clouds are located both in the disks and halos of nor- 
mal galaxies. In order to synthesize an absorption line profile, 
we run a line of sight through a model galaxy with a given 
orientation, impact parameter, and size (where the size of the 
galaxy is defined by its luminosity). We used the luminosity 
and impact parameter from Table g as the input for each galaxy 
model. The orientation of galaxy G5 was set to be edge-on but, 
when the other galaxies were used, their orientations were cho- 
sen randomly. Each galaxy model was chosen from 75D/25H 
Hybrid2 in Charlton & Churchill (1998), with 75% disk clouds 
and 25% halo clouds, by number. In order to simulate absorp- 
tion from either compact groups or galaxy pairs, the individual 
synthetic absorption profiles were shifted in velocity based on 
observed compact group or galaxy pair kinematics. We gener- 
ated 1000 realizations of synthetic spectra for each. 

4. 1 . Compact Groups 

The med ian velocity disper sion of compact groups, 
200 km s _1 (Hickson et al. 1992), results in profiles that have 
velocity spreads larger than in the z — 0.7443 system. A 
smaller velocity dispersion was required so that our models did 
not produce synthetic profiles with total velocity spreads incon- 
sistent with the data. Therefore, we chose the galaxy relative 
velocities from HGC 22, the c ompact group with th e minimum 



velocity spread (100 km s , Hickson et al. 1992). Thus, we 



have already restricted ourselves to a limited region of parame- 
ter space; this is not a "typical" compact group. 

Twenty-five randomly selected synthetic spectra from our 
three-galaxy compact group models (using G3, G4, and G5) 
are shown in the upper panel of Figure |[ Clearly, the velocity 
spreads of the clouds in the simulated spectra tend to be greater 
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than that in the z — 0.7443 system, even with the bias towards 
a small dispersion in galactic velocities. Also, the saturation 
levels in the simulated spectra tend to be higher, implying that 
our particular line of sight passes through less material than in 
an average model compact group. Finally, the regular struc- 
ture of pairs in absorption, separated by ~ 30 km s _1 , is not a 
characteristic feature of lines of sight through model compact 
groups. 

4.2. Galaxy Pairs 

In order to reduce the total velocity spread and saturation, we 
also considered models of galaxy pairs, assuming that either G3 
or G4 is at z — 0.7443 along with G5. Again, G5 was taken 
to have an edge-on orientation, and the other galaxy was taken 
to be randomly oriented. The velocity difference between the 
galaxy pair is tuned to ~ 50 km s _1 . In fact, this is a relatively 
typical value for field galaxy pairs. 

Twenty-five randomly selected galaxy pair synthetic spectra 
are displayed in the bottom panel of Figure ^. Compared to the 
compact group synthetic profiles, a larger percentage of the pair 
profiles have velocity spreads and saturation levels consistent 
with the observed Mg II profile. However, the regular absorp- 
tion line pairs we observe in the spectrum of MC 1331 + 170 
are only rarely produced, even in this simple galaxy pair model. 

4.3. Kinematic Model Results 

The z = 0.7443 absorber toward MC 1331 + 170 was se- 
lected for this study because of its uniquely regular properties. 
As such, this is a biased study. Nonetheless, we have seen that, 
in 1000 Monte-Carlo simulations of galaxy pairs and in 1000 
simulations of compact groups, we can only rarely produce a 
profile with such low levels of saturation and regular velocity 
splittings. However, it is difficult to quantify the shapes of the 
profiles in such a way that we can make direct statistical com- 
parisons. We remind the reader that we used a compact group 
with an unusually small velocity dispersion. Therefore, based 
on this and on a visual inspection, we conclude that the ob- 
served profile is unlikely to arise from either a compact group 
or galaxy pair environment. 

5. SUPERBUBBLE SIMULATIONS 

The spectrum of the z — 0.7443 absorption system appears 
to be kinematically consistent with a line of sight intersecting 
the oppositely expanding sides of th e "cold" shells of multiple 
superbubbles (Churchill et al. 1995). In order to test this hy- 
pothesis, we calculated what the column densities of Mg II and 
Fe II should be in the context of such a model. 

5.1. The Model 

The global properties of our model superbubbles are gov- 
erned by the similarity solut ions of Mac Low & McCray (1988) 
(also see Weaver et al. 1977 ). Models of the cold shell of the su- 



perbubble s were performed using Cloudy, version 94.00 (F er- 
land 1996), taking into account chemical and ionization con- 
ditions. The photoionizati on source consists of the ultraviolet 
background from quasars ( Haardt & Madau 1996[ ) and the O 
and B stars driving the expansion of the superbubble. The in- 
tensity of the Haardt-Madau background is set to a constant 
hydrogen-ionizing photon number density of 10~ 52 cm -3 (as 
appropriate for z ~ 1). The stellar contribution depends on 
the luminosity of the parent OB association and the radius of 
the superbubble. Based upon our models, we find that for lu- 
minosities greater than ~ 10 38 ergs s -1 , the stellar spectrum 



dominates at ages less than ~ 10 years. Most of the models, 
however, are dominated by the Haardt-Madau spectrum. 

We assume the ISM to have a constant initial atomic density 
of no cm' 3 . An OB association residing in this medium gen- 
erates a total luminosity of L^g (in 10 38 ergs cm -3 ) through 
a series of supernova explosions. If each supernova produces 
~ 10 51 ergs and the OB association lasts 50 million years, 
the time-averaged luminosity is related to the number of su- 
pernovae by N$n — 580 L^g. The atomic density inside of the 
swept up shell is n e (expressed in units of cm -3 ) and the metal- 
licity is Z, in solar units. The energy input from supernovae 
drives the b ubble radially outward into the ISM. The similar- 
ity solution (Mac Low & McCray 1988) gives the radius of the 
bubble, 



R = 0.27 {n / ^ 



r a2 t° 7 6 



kpc, 



(1) 



where tj is the age in units of 10 7 years. The shell expansion 
speed then follows as 



R = 16 (n /L 38 ) 



-0.2 



km s 



(2) 



which is constrained by the observed velocity splitting, Av, of 
the absorption lines according to 



R = 0.5Av (1 



2 /R 2 



-1/2 



(3) 



where z is the impact parameter (to the center of the bubble) 
of the line of sight. This equation is not very sensitive to im- 
pact parameter unless z/R ~ 1. From the observed profile, 
Av ~ 30 km s _1 ; thus, we set R = 15 km s _1 . 

If we choose L^s and tj, we can calculate R and N to t for 
input to the Cloudy model, where N to t is the total hydrogen 
column density. The radius is found by combining equations 
(@) and (§), giving 



R = 0.017 Rt 7 kpc. 
From equation ([!]), n is given by 

n Q = L 38 (l6if°' 4 AR) 5 cm 
Finally, N tot is 

N tot = 10 21 R n cm- 



-3 



(4) 



(5) 



(6) 



The luminosity, L^s (0.01,0.1,1,10), was chosen for OB as- 
sociations contai ning between 5 and 5, 000 stars (e.g. M cKee 
& Williams 1997; presolin et al. 1998[ ). The age, t 7 (0.5,1,2,4), 
was chosen to represent different points in the evolution of a 
bubb le expanding for 50 million years (M ac Low & McCray 
1988). The Cloudy models of the shell required that we also 
input logn e (—2.5 to 2, with intervals of 0.25). Cloudy was 
then run for every combination of the chosen parameters (L38, 
tj, n e , and Z). The metallicity ran ge (l og Z <J 0) was based 
upon observational constraints (see § 5.2). Each simulation was 
halted when the hydrogen column density (N to t) calculated in 
Equation || was reached. 
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5.2. Simulation Results 

For each Lz$ and Z, Cloudy grids were generated. In Fig- 
ure [I| we show a plot of output grids from Cloudy photoioniza- 
tion models for the \ogZ = — 1 and logZ = —0.5 cases. The 
thickness of the plane parallel supernova shell is plotted versus 
its Mg II column density. Each panel represents a different lu- 
minosity of the parent OB association (given by L^g). The solid 
lines represent the ages of the model superbubbles, correspond- 
ing (from left to right) to tj = 4, 2, 1, and 0.5. The dashed 
lines represent the hydrogen number density of the shell, which 
ranges from (top to bottom) logn e = —3 to \ogn e — 2 cm -3 
in logarithmic intervals of 0.25. We picked the regions of pa- 
rameter space which are in the observed Mg II column density 
range (see Table [0 ) and are consistent with observed superbub- 
bles ( Heiles 1979] ). These regions are enclosed by the dotted 
boxes. 

The models produce values for n e /no over a large range, the 
maximum being ~ 1000. We only plot models for which the 
shell is more dense than the ISM (n e > no, where no is found 
from Equation^). This will probably be the case unless there is 
an extre me ISM density gradient (M ac Low, McCray, & Nor- 
man 1989 ) or we are passing through a cloudy medium ( Silich 
et al. 1996). Real superbubbles have shell densities greater than 
four times the ISM dens ities, perhaps greater t han ten times if 
the shells are overstable (Ryu & Vishniac 1988). 

It may seem counter-intuitive that Mg II column density de- 
creases as time increases, but the successive tj grid lines do not 
represent an evolutionary sequence. The decrease is a conse- 
quence of holding shell speed constant for all models. In or- 
der to produce the same speed at a later time, the ISM density 
muchbe lower (providing less resistance). Thus, the different 
tj lines represent varying initial conditions (see Equation ||). If 
you have a smaller ISM density, then you will have a smaller 
amount of material swept up into the shell, resulting in a smaller 
Mg II column density. 

Increasing the luminosity of the clusters corresponds to a 
shift of the grid toward higher Mgll column densities and 
thicker shells. Increasing metallicity causes a shift in the grid 
towards lower Mgll column densities. We have set constraints 
on the value of L^s, so we can now constrain Z to fit ob- 
servational parameter space (the dotted box). The three left- 
hand panels of Figure H show log Z = — 1, which is consistent 
with the observed Mgll column densities for L^s = 0.01 or 
L38 = 0.1. A grid with an increased metallicity would require 
a lower L38 (as apparent in the three right-hand panels). If 
log Z > 0, the models require a luminosity that is too low for an 
OB association (L^s < 0.01). The l argest OB associations con - 
tain ~ 7000 supernova progenitors ( McKee & Williams 1997 ), 
corresponding to L38 ~ 10 (not shown in Figure Q). Even un- 
der these extreme conditions, we cannot place a useful lower 
limit on the metallicity in the shell. 

We can not obtain general lower and upper limits on tj 
(though tj < 5 is a physical limit in our models, see § 5.1). 



The ages are limited, however, for particular combinations of 
I/38 and Z. For example, the L38 = 0.1, logZ = — 1 grid in 
Figure |] has only 1 < tj < 4 lying within the allowed parame- 
ter space (dotted box). Presumably, the multiple bubbles giving 
rise to our spectrum could arise from the same burst of star for- 
mation and have similar ages and metallicities. For this to be 
true, a logarithmic spread of 1.5 is required in L38 in order to 
reproduce the entire range of Mg II column densities. 

As with t-j, we cannot obtain general lower and upper limits 



on the shell density, n e , but we can constrain it for particular 
combinations of L38 and Z. As can be seen in Figure |], in- 
creasing L38 raises the grid and brings higher shell densities 
into the allowed parameter space. For example, the L38 = 0.1, 
logZ = —1 grid has only —1.5 > logn e > 1 cm~ 3 lying 
within the allowed parameter space. 

6. CONCLUSION AND DISCUSSION 

The line pairs in the z = 0.7443 Mgll absorption system 
in the spectrum of MC 1331 + 170 are strongly suggestive of 
shells of expanding superbubbles. The identity of the host of the 
z = 0.7443 absorption system is not immediately obvious, but 
is likely to be the spiral galaxy (G5) seen in the //5T/WFPC2 
image (Figure ^). G3 and G4 are also candidates for being 
at z = 0.7443 and may form a compact group or pair with 
G5. However, kinematic models of galaxy pairs and compact 
groups seem to yield velocity spreads larger than the observed 
~ 200 km s _1 and rarely reproduce the symmetry seen in the 
MC 1331 + 170 system. This evidence seems to suggest that 
the absorption arises in a single galaxy, either from superbub- 
bles or random cloud distributions. Photoionization models, 
considering the conditions of the parent OB association and the 
surrounding medium, support the superbubble hypothesis with- 
out requiring fine-tuning of the specific conditions from which 
the superbubbles arise. 

The models of superbubbles do, in fact, produce the ob- 
served Mgll column densities and velocity splittings of the 
z = 0.7443 system, but there are still some discrepencies. 
The ~ 250 km s _1 spread of our profile is the most prob- 
lematic feature in the context of disk superbubbles. A typi- 
cal L* galaxy has a rotation velocity (Vq) of ~ 220 krns" 1 . 
If one assumes that the superbubbles follow the disk rotation, 
the maximum spread (V max ) over which they could occur is 
V m ax = Vb(l — D/R g ), where D is the impact parameter and 
R g is the maximum galactic radius at which superbubbles can 
occur. The parameter R g is likely to be approximately equal 
to the optical radiu s of the galaxy. In M101, superbubbles are 
seen out to 30 kpc (Kamphuis 1993). If we assume that R g can 
be as large as ~ 40 kpc, we get only V max ~ 100 km s _1 . For 
G5 at z = 0.7443, the impact parameter is 20h~ 1 1 kpc. 

There are several possible explanations for the large spread 
in the z = 0.7443 Mgll absorber: 1 ) A merger is inducing 



globular cluster forma tion in the halo (Ashman & Zepf 1992 



^epf & Ashman 1993) and the clusters are creating superbub- 
bles. 2) Superbub bles are forming within tidal debris fro m a 
recent interaction (Kniermann et al. 2000; Gallagher 2001). A 
line of sight through ti dal debris can cr eate velocity spreads as 
large as 1000 km s" 1 (Gallagher 2001 ). 3) Some of the super- 
bubbles are outliers in vel ocity space, possibly locate d within a 
companion dwarf galaxy (Mac Low & Ferrara 1999). 4) Some 
of the absorption lines are not from superbubbles and lie out- 
side of the disk. This absorption could arise in halo clouds or 
in dwarf galaxies. 

Womble, Junkkarinen, and Burbidge (1991) point out that 
large velocity spreads in edge-on galaxies are better explained 
by chimneys, worms, or superbubbles blowing out into the halo. 
These are processes known to occur at the end of a superbub- 
ble's lifetime when the shell accelerates through a density gra- 
dient into the halo. The shell then fragments due to Rayleigh- 
Taylor instabilities and the bubble begins expelling material out 
of the disk. This process is quite messy, however, and is not 
likely to give rise to such evenly split lines. 
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Our spectrum shows no more than six superbubbles within 
the line of si ght, each ~ 10 pc thick ( Mac Low, McCray, & 
Norman 1989); thus, there must be a substantial path length 
through the ISM in between the superbubbles. Why don't we 
see dominating absorption from the ISM between the super- 
bubbles? The simplest explanation is that the filling factor of 
the superbubbles is so large that the path length through the 
ISM makes a negligible contribution to t he absorption. O bser- 
vations of Hi holes in nearby galaxies (Kamphuis 1993) give 
maximum covering factors of ~ 0.3; however, they are only 
sensitive down to R ~ 0.5 kpc. Another possible explana- 
tion is that OB associations reside in regions of enhanced den- 
sity and the superbubbles sweep the densest material into their 
shell before reaching the sparse, intervening material. Finally, 
the parent OB associations themselves are at unusual velocities 



with respect to most of the ISM material. This separation in 
velocity space may arise if the bubbles are localized in a certain 
part of the disk or exist in a star-forming tidal tail nearby. The 
lines at ~ —100 km s _1 in Figure [j] show less symmetry than 
the others and could be arising from ISM clouds in the disk. 
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Ql 33 1 + 170 Field Analysis Scenarios 
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Fig. l.— The HIRES/Keck spectrum of the z=0.7443 system toward Q1331+170 showing transitions detected at > 5er aligned 
in rest frame velocity space. The solid lines represent a model from a simultaneous Voigt profile fit to the spectra, with the ticks 
marking individual components. 



BOND ET AL. 



9 




Fig. 2. — Two images, each with a different dynamic range, of the MC 1331 + 170 field combining the R (F702W) and I (F814W) 
filters. Galaxies are numbered G3 through G5, in order of distance from the quasar (at center). G2 is not labeled because of its 
proximity to the quasar and can only be resolved via PSF subtraction. The left-hand panel emphasizes G3 and G4 and the right-hand 
panel emphasizes G5 (for presentation purposes only). 
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Fig. 3. — A color-redshi f t diagram plotted of model spectral energy distributions representing a variety of morphological types 
(Bruzual & Chariot 1993; Kinney et al. 1996). The dotted vertical lines denote the redshifts of the absorption systems in the 
spectrum of Q1331+170 and the dashed horizontal lines denote the colors of the sources in the optical field. 
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Fig. 4.— Output grids from Cloudy photoionization models for the log Z = — 1 and log Z = —0.5 cases. Each plot of shell thickness 
vs. Mg II column density in a column presents a box that represents a different luminosity input from the parent OB association (given 
by L^%). From left to right, the solid curves represent the ages of the model superbubbles, corresponding to t-j = 4, 2, 1, and 0.5. The 
dashed curves represent hydrogen number densities in the shell which are (from top to bottom) log n e = —3 to 2 cm~ 3 in logarithmic 
intervals of 0.25. The boxes in the center of the plots represent expected parameter ranges based on observed superbubbles. 
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Fig. 5.— Spectra generated from 25 models of three-member compact groups (top panel) and 25 models of galaxy pairs (bottom 
panel). 



